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£\j ' Abstract 

■ If a body of dielectric material is coated by a plasmonic structure of negative dielectric 

D ' constant with nonzero loss parameter, then cloaking by anomalous localized resonance (CALR) 

^"*^ . may occur as the loss parameter tends to zero. The aim of this paper is to investigate this 

phenomenon in two and three dimensions when the coated structure is radial, and the core, 
shell and matrix are isotropic materials. In two dimensions, we show that if the real part of the 
permittivity of the shell is —1 (under the assumption that the permittivity of the background 
is 1), then CALR takes place. If it is different from —1, then CALR does not occur. In 
. three dimensions, we show that CALR does not occur. The analysis of this paper reveals that 

1 occurrence of CALR is determined by the eigenvalue distribution of the Neumann-Poincare- 

type operator associated with the structure. 

-5' 

1 Introduction 

If a body of dielectric material is coated by a plasmonic structure of negative dielectric constant 
(with nonzero loss parameter), then anomalous localized resonance may occur as the loss parameter 
tends to zero. This phenomena, first discovered by Nicorovici, McPhcdran and Milton [36] (see 
also [33]), is responsible for the subwavelength focussing properties of super lenses [35], and also 
occurs in magnetoelectric and thermoelectric systems [33J . The fields blow-up in a localized region, 
l/"") , which moves as the position of the source is moved, which is why it is termed anomalous localized 

resonance. Remarkably, as found by Milton and Nicorovici [30] the localized resonant fields created 
by a source can act back on the source and cloak it. This invisibility cloaking has attracted much 

(N ' attention [30] 137] E |3T]|3S]|31 [571 QHI3J [2H3U3]. 

To state the problem and results in a precise way, let f2 be a bounded domain in M. d , d = 2, 3, 
and D be a domain whose closure is contained in f2. For a given loss parameter S > 0, the 
permittivity distribution in M. d is given by 



in R d \ H, 

e s ={e s + iS inn\D, (1.1) 
in D, 
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where —e s and e c are positive. We may consider the configuration as a core with permittivity e c 
coated by the shell Q \ D with permittivity e s + i8. For a given function / compactly supported 
in M. d \ satisfying 



fdx = (1.2) 

(which is required by conservation of charge), we consider the following dielectric problem: 

V-egVVg = f inR d , (1.3) 

with the decay condition Vg{x) — > as |x| — > oo. The problem of cloaking by anomalous localized 
resonance (CALR) can be formulated as the problem of identifying the sources / such that first 

E s :=[ S\\7V s \ 2 dx -> oo as 6 ->• 0, (1.4) 
Jn\D 

and second, Vg/ ' \[Eg goes to zero outside some radius a, as 6 — > 0: 

\Vg{x) / y/E~g\ -> as 5^0 when|x|>a. (1.5) 

Physically the quantity Eg is proportional to the electromagnetic power dissipated into heat by the 
time harmonic electrical field averaged over time. Using integration by parts we have the identity 

Eg = %[ (e 5 VV s )-VV s dx=-$ [ fV s dx (1.6) 

which equates the power dissipated into heat with the electromagnetic power produced by the 
source, where Vg is the complex conjugate of Vg. Hence (|1.4|) implies an infinite amount of energy 
dissipated per unit time in the limit S — > which is unphysical. If we rescale the source / by 
a factor of l/y/Eg then the source will produce the same power independent of S and the new 
associated potential Vg/y/Eg will, by (|1.5[) . approach zero outside the radius a: cloaking due to 
anomalous localized resonance (CALR) occurs. 

In the recent paper [2] the authors develop a spectral approach to analyze the CALR phe- 
nomenon. In particular, they show that if D and Q are concentric disks in R 2 and e c = — e s = 1, 
then there is a critical radius r* such that for any source / supported outside r* CALR does not 
occur, and for sources / satisfying a mild condition CALR takes place. The critical radius r* is 
given by 

r* = v^lM, (1-7) 

where r e and are the radii of H and D, respectively. It is worth mentioning that these results 
were extended in to the case when the core D is not radial by a different method based on a 
variational approach. 

The purpose of this paper is to extend some of the results in [2] in two directions. We consider 
the case when e c and — e s are not both 1 and we consider CALR in three dimensions. The results 
of this paper arc threefold: Let f2 and D be concentric disks or balls in M. d of radii r e and r,, 
respectively. Then, the following results hold: 

• If d = 2 and e s = —1, then CALR occurs. When e c = 1 the critical radius r* is given by 
(|1.7p and when e c ^ 1 the critical radius is 

r 2 

r. = f. (1.8) 

That is, for almost any source / supported inside r» CALR occurs and for any source / 
supported outside r* CALR does not occur. When e c 7^ 1 the cloaking radius r 2 /r^ matches 
that found in |30j for a single dipolar source (see figure 5 there and accompanying text). 
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• If e s — 1, then CALR does not occur. 

• If d = 3, then CALR does not occur whatever e s and e c are. 

We emphasize that the result on non-occurrence of CALR in three dimensions holds only when 
the dielectric constant e s is constant. In the recent work [3] we show that CALR does occur in 
three dimensions if we use a shell with non-constant (anisotropic) dielectric constant. 

It turns out that the occurrence of CALR depends on the distribution of eigenvalues of the 
Neumann-Poincare (NP) operator associated with the structure (see the next section for the defi- 
nition of the NP operator). The NP operator is compact with its eigenvalues accumulating towards 
0. It is proved in [2] that in two dimensions the NP operator associated with the circular structure 
has the eigenvalues ±p" for n = 1,2,..., where p = Ti/r e . We show that in three dimensions the 
NP operator associated with the spherical structure has the eigenvalues 

± 2 (2n + l) A/l + 4 " ( " +1)p2n+1 ' n = ( L9 ) 

The exponential convergence of the eigenvalues in two dimensions is responsible for the occurrence 
of CALR and the slow convergence (at the rate 1/n) in three dimensions is responsible for the 
non-occurrence. 



2 Layer potential formulation 

Let G be the fundamental solution to the Laplacian in R d which is given by 



G(x) 



^mM, d = 2, 



Let Ti := dD and T e := dfl. For F = T.; or T e , we denote the single layer potential of a function 
<p € L 2 (T) as SrM, where 

Sr[(p](x) := J^G(x-y)cp(y)da(y), 

We also define the boundary integral operator A^r on L 2 (T) by 

£r[<p](x) := [ dG ^. X , V \ {y)d<T{y), x G F, 
Jr ou(y) 

and let /Cp be the L 2 -adjoint of ICr- Hence, the operator /Cp is given by 

dG(x-y) ^ t2 , 



r 



The operators Kt and /Cp are sometimes called Neumann-Poincare operators. These operators are 
compact in L 2 (T) if F is C 1,a for some a > 0. 

The following notation will be used throughout this paper. For a function u defined on R d \ T, 
we denote 



and 

du 
dv 



u\±(x) := lim u(x±tv(x)), iGT, 
t->o+ 



(x) := lim (Vu(x ± tv(x)), v(x)) , x £ T, 
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if the limits exist. Here and throughout this paper, ( , ) denotes the scalar product on R d . 

The following jump formula relates the traces of the normal derivative of the single layer 
potential to the operator /C r . We have 



^Srl^Jx) = (±±i + Kz\[p](x), xer. 

Here, v is the outward unit normal vector field to T. See, for example, [51 112|. 
Let F be the Newtonian potential of /, i.e., 



(2.1) 



F{x) 



G(x-y)f(y)dy, xeR d . 



(2.2) 



Then F satisfies AF = f in K d , and the solution Vs to (|1.3|) may be represented as 

V s (x) = F(x)+Sr t [<Pi](x)+Sr t [<p e ](x) 



(2.3) 



for some functions (pi <E L^iTi) and tp e £ L^{T e ) (Lq is the collection of all square integrable 
functions with zero mean- value). The transmission conditions along the interfaces T e and I\ 
satisfied by Vs read 



(e s + 15)-^— 
dv 



dVs_ 
dv 



dV s 
dv 

[e s + td) — 
dv 



on Ti 

on r„ 



Hence the pair of potentials (tpi,(p e ) is the solution to the following system of integral equations: 



(e s + id) 



dS Ti [*Pi 



+ (e s 



id) 



dS Tc [t£e 



dF 

(— e s + e c - id)—— on T, 
ov. 



, , , ■ c\ 9S Ti [cpi] dS Te [<Pe] 
(-1 + e s + id) 



Oi 



dv K 



+ (e s + id) — 5 

+ ov e 



= (1 - e s - id) 



dF 

dv e . 



on r R . 



Note that we have used the notation Vi and v e to indicate the outward normal on Ti and T e , 
respectively. Using the jump formula (|2.1[) for the normal derivative of the single layer potentials, 
the above equations can be rewritten as 



zfl-K 
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r, 



Sr 



. dv, 



■S 



r, 







- dF - 








dvi 








dF 






- dv e _ 



(2.4) 



on Uq = Ll(Ti) x Ll(T e ), where we set 



e c + e s + id 
2(e c - e s - id) ' 



l + e s + id 
2(1 -e a - id)' 



(2.5) 



Let 7-L = L 2 (Ti) x L 2 (T e ) and let the Neumann- Poincare- type operator K* : % — > % be defined 



by 



d 

■Apr. — - — <Sr c 



,dv e 



dvi 

IK" * 



(2.6) 
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and let 



$ := 



(fie 



dF - 

dvi 
dF 
dv K . 



Then, (|2.4[) can be rewritten in the form 



where I s is given by 



(I A '+K*)$ = .g, 



zfl 
zU 



(2.7) 

(2.8) 
(2.9) 



3 Eigenvalues of the NP operator 

It is proved in [2] that for arbitrary-shaped domains il and D the spectrum of the NP operator K* 
lies in [—1/2, 1/2], and if fi and D are concentric disks, the eigenvalues of K* on Ho are ±p n /2, 
n = 1, 2, . . .. In this section we compute the eigenvalues K* on H when Q, and D are concentric 
disks or balls. 

3.1 Two dimensions 

Let r = {|x| = ro} in two dimensions. It is known that for each integer n 



S r [e m °](x) 



ro 


- 


2\n\ 




ro 




2\n\ 





An6 



if |x| = r < r , 
if |a;| = r > tq. 



Moreover, 
and 



/Cr[l] = i 



(3.1) 

(3.2) 
(3.3) 



In other words, JCr is a rank 1 operator whose only non-zero eigenvalue is 1/2. 

Using (|3.2p . it is proved that eigenvalues of K* on Ho are ±p 2 /2 (see [5]). We now show that 
±1/2 are also eigenvalues of K* on Ho- These eigenvalues are of interest in relation to estimation 
of stress concentration [3]. Using (|3.3[) we have 



Sr[l](aO 




if 



r < r , 



and hence 



It then follows that 



dr 



5r[l](a;) = <^ l 



if \x\ = r > r , 

if |x| = r < r , 
if |x| = r > ro- 








r-i o- 




a 


& 




i i 




b 



(3.4) 



(3.5) 



(3.6) 



where a and 6 are constants. So ±1/2 are eigenvalues of K* 
We summarize our findings in the following proposition. 
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Proposition 3.1 The eigenvalues o/K* defined on concentric circles in two dimensions are 



and corresponding cigenf unctions are 



1 




"0" 




" ±in8 - 




" e ±in6 ' 


_ J_ 

T e _ 




1 


) 


pe ±in6 


1 


_ pe ±in0 



n = 1,2, 



(3.7) 



(3.8) 



3.2 Three dimensions 

Let Y™(x) (to = — n, —n + 1, . . . , 0, 1, . . . , n) be the orthonormal spherical harmonics of degree n. 
Here x = A. Then |a;| n l^"(x) is harmonic in M 3 . 

Lemma 3.2 Let T = {\x\= ro} in three dimensions. We have for n = 0, 1, . . . 

1 



2(2n+ 1 

Proof. It is proved in [TSJ Lemma 2.3] that 



Y^(x), M =r "o, m = -n,...,n. 



(3.9) 



2r 

for any function (p <G L 2 (T). So it follows from (|2.1[) that 

(9 11 

— 5 r M|_(x) + —S r [<p]{x) = --<p(x), \x\ = ro- 
ar 1 2ro 2 



(3.10) 



(3.11) 



Let f{x) = Y™{x). Then since Sr[Y™](x) and |x| n y^™(i) are harmonic functions in {\x\ < ro}, 
we have 

1 r n 



Sr[Y™](x) = 



2n 



-—Y™{x), for \x\ = r < r , 
1 r n 



(3.12) 



and (JSU) follows from ([3~T0)l . □ 
Lemma l3.2l savs that the eigenvalues of /C r on L 2 (T) when T is a sphere are 2 (2n+i) » n = 0, 1, . . ., 

and their multiplicities are 2n + 1. 
By (j3"7T2"]) . we have 



2n + 1 V r 



ra-1 



Y™(x), \x\= n 



(3.13) 



Similarly, we have 



and hence 



n+2 



We now have for constants a and b 

b^- T p n - 1 )Y: n 







bY m 

_ n 





2(2ra+l) T u 2n+l P 

n+1 „n+2 I 6 
"2ri+l^ ^ 2(2n+l) / " 



2(2?i+l) 2)1+1 P 
1 



n+1 rt+2 
L2n+1^ 



(3.14) 



(3.15) 



Thus we have the following result. 
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Proposition 3.3 The eigenvalues of K* defined on two concentric spheres are 

1 



± olo ^.^ 1 + 471(71 + 1)^+1, n = ,l,..., (3.16) 
2(2n + 1) 



and corresponding eigenfunctions are 

~t /l i i 1 ^ „2n+l 1\\rml 

771 = —77, . . . , 77, (3.17) 



(v/i + 477(77 + i)p 2 "+ 1 - i)y n * 
2(71 + i)p™+ 2 y I ™ 



(-^1 + 477(77 + i)p2»+ 1 - i)y n r 
2(71 + i)p™+ 2 y™ 



respectively. 



It is quite interesting to observe that if we let | = Ao > Ai > . . . be the eigenvalues of K-r for a 
disk or a sphere enumerated according to their multiplicities, then the eigenvalues \x n of K* satisfy 

fi n = ±\ n + 0{p n ). (3.18) 

4 Anomalous localized resonance in two dimensions 

In this section we consider the CALR when the domains f2 and D are concentric disks. We first 
observe that zf and z% converges to non-zero numbers as 5 tends to if e c 7^ — e s 7^ 1. So, in this 
case CALR does not occur regardless of the location of the source. Furthermore, if e c = e s = 1, 
a thorough study was done in [2J. It is proved in [5] that if the source / is supported inside the 
critical radius r* = y/rf/ri, then the weak CALR occurs, namely, 

limsupE'a =00. (4.1) 

Moreover, if F is the Newtonian potential of / and the Fourier coefficients <?" of —■§£- satisfies the 
following gap property: 



[GP] There exists a sequence {nk} with \m\ < \n 2 \ < ■ ■ ■ such that 

\ n n k 12 

lim pK+d-KI Igg = 00; 

fc-s-oo |71fe|pl nfc 



then CALR occurs, namely 



limi^oo, (4.2) 

5— >0 



and Vs/y/Es goes to zero outside the radius yrf/r; • 

The remaining two cases are when e c 7^ — e s = 1 and e c = — e s =/= 1. In these cases, we have the 
following theorem. 

Theorem 4.1 (%) If e c = — e s 7^ 1, then CALR does not occur, i.e., 

E 5 <C (4.3) 

for some C > 0. fVFe note, however, that there will be CALR for appropriately placed sources 
inside the core, as can be seen from the fact that the equations are invariant under conformal 
transformations, and in particular under the inverse transformation 1/z where z = X\ + ix 2 , 
which in effect interchanges the roles of the matrix and core.) 

(ii) If e c 7^ — e s = 1, then weak CALR occurs and the critical radius is 7% = t\t^ , i.e., if the 
source function is supported inside r* ( and its Newtonian potential does not extend harmon- 
ically to R 2 J, then 

lim sup = 00, (4.4) 
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and there exists a constant C such that 

\V s (x)\<C (4.5) 

for all x with \x\ > r\jr\. 

(Hi) In addition to the assumptions of (ii), the Fourier coefficients g™ of — satisfies the fol- 
lowing gap property: 

[GP2] There exists a sequence {n^} with \n\\ < \n2\ < ■ ■ ■ such that 

\ n n k 1 2 

lim p 2(|n fc + 1 |-|n fc |) \9e ^ ^ 

lim E s = oo, (4.6) 

i—5-0 

and Vg/^/Es goes to zero outside the radius r^/rf, as implied by l[4.5\ ). 

Before proving Theorem 14.11 we make a remark on the Gap Properties [GP] and [GP2]. One 
can easily see that [GP] is weaker than [GP2], namely, if [GP] holds, so does [GP2]. 

The rest of this section is devoted to the proof of Theorem 14.11 As was proved in [5] , we have 

— SrAe m6 ](x) = -p M+1 e me , 

^ l Sr e [e me ](x) = -ipl"l- 1 e me . 

Using these identities, one can see that if g denned by (|2.7p has the Fourier series expansion 



then the CALR occurs, i.e., 



■* = £ 



e ine 



then the integral equations (|2.8[) are equivalent to 

A P H + 

Z e <Pe-\ —ft = 9e 

for every |n| > 1. It is readily seen that the solution $ = (c/5,-, ip e ) to (|4.7|) is given by 



<Pi = 2 2^ —TT-x I^^ e 



Azfzf - p 2 ! r 



= 2£ 



2zf ge " - ^"l+^r 
4z/a£ - p 2 l"l 



If the source is located outside the structure, i.e., / is supported in |a;| > r e , then the Newtonian 
potential of /, F, is harmonic in |a;| < r e and 

F(x)=c-J2— i tz I r ]nl e ine , \x\<r e . (4.8) 



n#0 



\n\r e 



Thus we have 

ff? = -fl?P |n|_1 - (4-9) 



So we have 



oV (2zf + l)pH-i g n 
(2^ + p 2 l"l)(7™ 



^ = 2 E 4zV- P " ; - 



Therefore, from (|3.1[) we find that 
and 



o/„2[n| A 2|n| „ 

2(r, z e r e zj ff ^ re<r = | x | 



i i \n\-X, 



Inlri™ 1 " 1 ^ 2 !"! -4zfzf)^ |n| 



r, < r = \x\ < r e 



^oH^ hl (p 2 l"'-4z^) 



ffe "rl"le me , < r = \x\ < r e 



We obtain the following lemma. 
Lemma 4.2 There exists 5q such that 



E 5 « ^ 



E 
E 



ii 2 



n|(,5 2 +p 4 l"l) : 



i/e c 7^ e s = 1, 



n |(<5 2 +/H)' 



ife c = e s ^ 1, 



uniformly in 5 < Sq. 

Proof. Using (|4~T2")) . and (|4TT^| . one can see that 



^(x) = c + r e 

n#0 



„ 2 I' 



■ |n| 



2z?r |n| 



(22* + l)5?e 



n „in8 



\n\rl nl {4zfz 5 e - pW) 



We check that 



2|n| 

Y 

r |n| 



2z ? r N e ™e 



2|» 



2|«| 



2zjVl r 



Then straightforward computations yield that 



/ 6\VV S \ 2 *J2 S 



2zt + 1 



(4|zf| 2 +p 2 H)M 



'BAB, „ #0 

If (5 is sufficiently small, then one can also easily show that 

\4zfz s e -p 2 W\K5 + p 2 W. 



n\ ■ 



Therefore we get (|4.14[) and the proof is complete. 
First, if e c = — e s =/= 1, then 



n |(J2 + p4|n|) - 2^ 2 |n| - 2 



<9f 



dv K 



L2 ( r e ) 



<C||/IU»(R» 



(4.10) 



(4.11) 

(4.12) 
(4.13) 



(4.14) 



□ 
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Suppose that e c 7^ — e s = 1, and let 



If \n\ < Ng, then 8 < /9 2 ' n ', and hence 



Ns = (4.15) 
21ogp 



y W > y W > v «1 (4i6) 



If the following holds 



lo n l 2 



then one can show as in [2] that there is a sequence {|fifc|} such that 



lim E ,„ fc = 00. (4.18) 

Suppose that the source function / is supported inside the critical radius r* = r^r" 1 (and 
outside r e ). Then its Newtonian potential F cannot be extended harmonically in \x\ < r* in 
general. So, if F is given by 

F = c- J2a n r^e me , r<r e + e (4.19) 

ra#0 

for some e > 0, then the radius of convergence of the series is less than r*. Thus we have 

limsup|a„| 2 r 2H = 00. (4.20) 



Since = \n\a n r^ \ (l4~T7ll holds. 
By (|4.11J) . we know 



|n| — >oo 



\V 5 \<\F\ + Cy - ^ <\F\+Cy 'J—^l<C (4.21) 

I "I — I I § + p2\n\ r |n| - I I 2|ra| r |n| — > 

n#0 r n=£0 ' i 

if r > /rf . Thus (ii) is proved. 

We now prove (iii). We emphasize that [GP2] implies (|4.17l) . but the converse may not be true. 
On the other hand [GP2] holds if 

lo n l 2 

lim Jpl = 00 . (4.22) 

So we may regard the condition [GP2] something between (|4.17[) and (|4.22l) . 

Suppose that [GP2] holds. If we take 6 = p 2a and let k(a) be the number such that 

K(a)| <a< K( Q ) + i|, 

then 

V J^l__- 2a J9ef_ > 2(|n fc(Q)+1 |-K M |) Iffe^T > f 4 95) 

0^|n|<JVj 1 0#|ri|<a 1 IP l'^(a)IF 

as a — ► 00. Combined with Lemma 14.21 and f|4. 16[) . it gives us (iii). 
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5 Non-occurrence of CALR in 3D 



In this section wc show that CALR does not occur in a radially symmetric three dimensional coated 
sphere structure when the core, matrix and shell are isotropic. We have the following theorem. 

Theorem 5.1 Suppose that F e and I\ are concentric spheres. For any e c and e s , there is a 
constant C independent of 5 such that if Vg is the solution to il.S\) . then 

S\VV S \ 2 < C||/||| 2(R3) . (5.1) 

n\D 

Proof. Suppose that -g^F has the Fourier series expansion 

t~. oc n 

^=-EE «■ (5-2) 

e n— m=—n 

Then one can show as in (14.90 that 



oo n 



i f = -EE fmnP^Y-. (5.3) 



n— m— — n 



By solving the integral equation (|2.4p using (|3.15[) . wc obtain 

oo n / 1 \ 

- E E A » V"" 1 K + « .9™„^ n , (5-4) 



oo n 

<A = - 

n— m— — n 

co n 



Y^ Y^ A -1 ( ^ 71 + 1 2n+l \ „e ym r\ 

^ = EL A »P ( * ~ 2(2^TT) + 2(2n~+T) P I 9mnYn ' (5 ' 5) 



n— m— — n 

where 



V 2(2n + l)J^ e + 2(2n+l)y' + (2n+l) 2P ' 
Suppose for simplicity that e c = — e s = 1, so that and given by (|2.5[) simplify to 



2(2 -i5)' 

Then one can see that if d is sufficiently small, then 

|A„| ^6 2 + n- 2 . 

So we have 



oo n 



s\\<Pi\\i HTi) <cj2 E 



V , e ,2 



^ (<5 2 + n- 2 ) 2 ' ffmTl1 

n— m— — n v J 
oo n 

H n 2ni„e |2 



<^EE ^P 2n \aU 

ri—0 m— — n 
oo 

< |^_| 2 < C||/||| 2(R3) . 



n=0 



and 



oo n X 2n 

Sy e \\h { r e) <Cj2 E s^- 2 \9 e m n\ 2 <C\\f\\h m . 



n—O m— — n 
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Therefore we have 

/ 6\VV S \ 2 = [ S\WF\ 2 + f 5\W(SrA^]+SrA^e})\ 2 
Jn\D Jn\D Jn\D 

<f S\WF\ 2 + f 5\W(SrA^}+SrA^e})\ 2 
Jn\D Jn\D 

< I 5\VF\ 2 + <y(||^||i a(r<) + ||Ve|li a( r.)) < C||/ll!=(H3)- 
Jn\D 

If e s ^ — 1 or/and e c =/= 1, then the same argument can be applied to arrive at (|5.1|) . This completes 
the proof. □ 
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